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Abstract. If a (cusped) surface S admits an ideal triangulation 
T with no shears, we show an efficient algorithm to give S as 
a quotient of hypebolic plane H 2 by a subgroup of PSL(2, Z). 
The algorithm runs in time O(nlogn), where n is the number 
of triangles in the triangulation T . The algorithm generalizes to 
producing fundamental groups of general surfaces and geometric 
manifolds of higher dimension. 



1. Introduction 

Let S be a cusped hyperbolic surface admitting an ideal triangula- 
tion with the following property: 

Property 1. For any pair of adjacent ideal triangles ABC and ABD, 
the cross ratio of the four points A, B, C, D equals 1, where the cross 
ratio is defined as: 

\A B C D] - (A - C)(B " D) 
[A ' B ' QD] - (B-C)(A-DY 

and we have implicitely identified the hyperbolic plane with the 
upper half plane HcC. 

This property has a number of equivalent formulations. One is 
geometric: 

Property 2. We can choose a family of horocycles h\,. . . , h n , where 
each hi is center on the z'-th cusp of of S, and hi is tangent to hi whenever 
Ci is adjacent to c ; in T. 
Another is algebraic: 
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Property 3. S is the quotient of H by a subgroup T of PSL(2, Z). 

In this note we will prove that the three properties are equivalent, 
and also give an algorithm to express S = r\H. The algorithm runs 
in time bounded by 0(n log n), and produces a set G of independent 
matrix generators for T. Since S is cusped, T is a free group, and so 
this is a complete description of T. In fact, we construct the generators 
as words in the two linear fractional transformations L and R, where: 

L(z) = z + 1, R(z) = — 

z — 1 

and the matrix generators are obtained by multiplying the words out. 
The algorithm consists of a number of steps: 

Step 1. Construct the Poincare dual T* of the triangulation. This will 
have a vertex for each face of T and a face for each vertex of T . This 
is an oriented complex, and thus we can cyclically order the edges at 
each vertex. 

Step 2. Construct a maximal spanning tree M of the 1-skeleton T* of 
T*. The edges of 7~* fall into two types. The edges of the first type, 
are the edges of M, the edges of the second type are not. 

Step 3. Split each edge of the second type. By "split", I mean that we 
replace the edge AB by a pair of edges AC\, BC2. We will henceforth 
refer to C\ and C2 as twins. 

After we split all the edges of the second type in the graph 7~*, 
we obtain a graph B, which is a tree where every non-leaf node has 
degree three. In addition, every leaf node is annotated with a cyclic 
ordering of the three edges. We are ready for: 

Step 4. Construct the shortest path from each leaf node to its twin. 
This path will look like C\V\ . . . VkCj. At each vertex Vk we have a fork 
in the road, and we annotate vy_ with an L or an R depending on 
whether we go left or right at the fork. 

Now we are done: each path from C\ to C2 gives a generator of the 
fundamental group of S, if we replace L and R by the linear fractional 
transformations with the same names (this should be done as we 
are constructing the paths, doing it after will bring us back to 0(n 2 ) 
running time). 

The correctness of the algorithm above follows immediately from 
the Poincare Polygon Theorem (see, eg, IITlO. 

1.1. Crossratios are all 1 if and only if there is a horodisk packing. 

. This follows from the observation that there is a unique horodisk 
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Figure 1. The modular tessellation (in the Klein model) 



packing of an ideal triangle. Indeed, if represent the ideal trian- 
gle ABC as one whose vertices are the three roots of unity in the 
Poincare disk model, the symmetric arrangement of horodisks obvi- 
ously works. Let the points of tangency of the horocycles (which are 
on the sides of ABC) be p A B, Pac, and Pbc- Now, suppose that there is 
another arrangements, with points of tangency q AB , CfAC, a BC, and let 
d(pAB, ^ab) = ?ab, and similarly for the other two sides. Suppose p AB 
lies between A and cj AB . Then the same is true of p AC and q A c- But 
the last two assertions would imply the the three new horocycles are 
not actually tangent along BC. To show the result we now need the 
following easy lemma: 

Lemma 1. Let T\ = ABC and T 2 = ABD be two adjacent ideal triangles, 
and let y be a horocycle centered on A. Let y x = y n T\, and yi = Y n T 2 . 
Then 

P\ = exp([A,B,C,D]). 

tyi I 

Proof. Let C = -1, A = oo, B = 0, D = z, and compute. □ 

1.2. All crossratios are 1 implies that the surface is a quotient of 
the upper halfplane by a subgroup of the modular group. This is 
not hard.to see, especially if one looks at the modular figure: each 
adjacent pair of colored and white triangles forms a fundamental 
domain for the action of PSL(2, Z) on the hyperbolic plane. Since 
the crossratios are all 1, the baricentric subdivisions of pairs of adja- 
cent ideal triangles agree, and so we see that our surface covers the 
modular orbifold. 
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2. Complexity 

2.1. Constructing the oriented dual (Step 1). The complexity of Step 
1 (constructing the oriented dual) depends on how one is given the 
triangulation. The most natural way is for it to be given as a rotation 
system, which is simply the graph with a cyclic ordering of the edges 
at every vertex. It is easy to see that in this case the dual graph can 
be constructed in time linear in the number of edges (the algorithm 
is simple: maintain a list of edges. Each edge is marked by or 
1. Initially all the edges have label 0. We pick the first edge e, and 
construct a list of edges obtained by always picking the edge which 
precedes e in the cyclic order. A closed cycle gives us a face (already 
equipped with the cyclic ordering of boundary edges). Every time 
an edge is seen we increase the label by 1. If the label is 2, we delete 
the edge from the list. Since each edge is seen at most twice, and we 
do constant work per edge, the algorithm is linear). 

The spanning tree (Step 2) can be done in time linear in the number 
of edges (see, eg, 0), and Step 3 can obviously be done in time linear 
in the number of vertices. This leaves us with Step 4, which we 
analyze below. 

2.2. Constructing the generators. At this point we have a tree (with 
every interior node of degree 3) M and a collection of pairs of leaves 
of M, and we need to construct paths between the two vertices in 
each pair. Since a shortest path between two vertices of a tree can 
be constructed in time 0(V(M)), (see Q) and the number of pairs 
is half the number of all leaf nodes (so 0(V(M)) as well), this gives 
an 0(V 2 (M)) algorithm for computing all the generators. We can do 
better, however, by first showing the following: 

Lemma 2. Let G be a tree with every non-leaf node having degree three. 
For every non-leaf node v, removing v separates G into three subgraphs 
G max (v),G med (v),G min (v),.with\V(G max (v))\ > |^(G med (z;))| > \V(G mhx (v))l 
Let v be the vertex which minimizes \V(G max (v))\. Then 

(1) V(G) 3 ~ 1 \ < V(G max (v))\ < ^V{G) + 1. 

Proof. Denote the three orders by N max , N med , N m m- The first inequality 
is true at any vertex, since 

(2) N max + N med + N med + 1 = W(G)\. 

To show the second inequality, let V\ be the vertex in G max (S) adjacent 
to v. At V\ the orders of the three components into which V\ separates 
G are N lf N 2f N 3 , where N 3 = N min + N med + 1, and N 1 +N 2 + l= N max . 
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Figure 2. Cut point 



Assume that N\ > N 2 - We then have two possibilities. The first is 
that N 3 > Ni. In that case, N med + N min + 1 > (N max - l)/2. Adding 
N max to both sides, we get the second side of the inequality!!] The 
second possibility is that N\ > N 3 . However, since Ni < N max , this 
contradicts the defining property of v. □ 

The next result we will need is the result of 0: Given a rooted 
tree T with a positive weight associated with every node, there is a 
linear time algorithm to partition the tree into a minimal collection 
of subtrees such that the weight of no subtree exceeds k. In our 
application, all the weights are equal to 1, and k = |V(T) + 2. By 
Lemma EJthe tree will be broken up into exactly two components. 

The algorithm is then simple: We partiion our tree into two rooted 
subtrees (the roots will be the two endpoints of the edge we delete to 
partition), For each of the two pieces we compute the pair distances, 
and all the distances from the leaves to the root (recursively), then 
use the distances we had computed to compute all the distances in 
the original tree. To make the first step the same as all the others, we 
pick an arbitrary leaf and call it the root. 

It is clear that at each step we have the following recurrence in- 
equality for the running time: 

T < cV(G) + T(Vi) + T(V 2 ), 

where T x + T 2 = V(G) - 1, and max V x , V 2 < §V(G) + 2. 
This clearly implies an 0(n log n) running time. 
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3. Extensions 

The algorithm described above easily extends to other cases. The 
simplest extension is where the ideally triangulated surface does not 
have all cross ratios equal to one. In that case, we simply replace the 
linear fractional transformations L and R by the appropriate conju- 
gates of the transformation L a/b (z) = az + b. 

For a surface (or a higher dimensional manifold) equipped with a 
triangulation by finite triangles, we simply develop the fundamental 
domain (as given by the spanning tree) into the model space, and 
then use the side-pairing information to produce the generators of 
the fundamental group. It should be noted that if a surface is finely 
triangulated, the generating set will contain many instances of the 
identity element, but this is obviously not a serious problem. 

It should be noted that Step 4 of our algorithm can be replaced by 
using the spanning tree information to embed the triangulation in H 2 , 
and then computing the relevant isometries. This would, however, 
necessitate running a version of the continued fraction algorithm for 
each side pairing, and thus will be much less efficient (0(n 2 ) vs our 
0(n logn).) 
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